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Abstract 

Transport events in turbulent tokamak plasmas often exhibit non-local or non-diffusive action 
at a distance features that so far have eluded a conclusive theoretical description. In this paper 
a theory of non-local transport is investigated through a Fokker-Planck equation with fractional 
velocity derivatives. A dispersion relation for density gradient driven linear drift modes is derived 
including the effects of the fractional velocity derivative in the Fokker-Planck equation. It is 
found that a small deviation (a few percent) from the Maxwellian distribution function alters the 
dispersion relation such that the growth rates are substantially increased and thereby may cause 
enhanced levels of transport. 

PACS numbers: 52.25.Dg,52.30.Gz,52.35.Kt 
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INTRODUCTION 



Understanding anomalous transport in magnetically confined plasmas is an outstanding 
issue in controlled fusion research. A satisfactorily understanding of the non-local features 
as well as the non-Gaussian probability distribution functions (PDFs) found in experimen- 
tal measurements of particle and heat fluxes is still lacking. In particular, experimental 
observations of the edge turbulence in the fusion devices [V\ show that in the Scrape of 
Layer (SOL) the plasma fluctuations are characterized by non-Gaussian PDFs. It has been 
recognized that the nature of the cross-field transport through the SOL is dominated by 

turbulence with a significant ballistic or non-local component where a diffusive description 

f] fl 
is improper [21]. Moreover, the scaling of the confinement time r oc with a < 2 \m is 

typical in low-confinement mode discharges, instead of the diffusion induced result t (x L"^, 

where L is the system size. There is a considerable amount of experimental evidence 0-0] 

and recent numerical gyrokinetic 0, [lH and fluid 12 1 simulations that plasma turbulence 



in tokamaks is highly non-local. 

In addition, intermittent turbulence is characterized by patchy spatial structure that 
is bursty in time. The PDFs of these intermittent events shows unimodal structure with 
"elevated" tails that deviates from a Gaussian prediction. The understanding of these events 
are at best limited lisll 



7| . Moreover, the high possibility of confinement degradation by 
intermittency strongly calls for a predictive theory. 

A prominent candidate for explaining the suggestive non-local features of plasma turbu- 
lence is the inclusion of a fractional velocity derivative in the Fokker-Planck (FP) equation 
leading to an inherently non-local description as well as giving rise to non-Gaussian PDFs of 
e.g. densities and heat flux. The non-locality is introduced through the integral description 
of the fractional derivative 18142 Ol| and the non-Maxwellian distribution function drives the 
observed PDFs of densities and heat flux far from Gaussian. 

The aim of this study is to elucidate the effects of a non-Maxwellian distribution function 
induced by the fractional velocity derivative in the Fokker-Planck equation. Some previous 
papers on plasma transport have used models including a fractional derivative where the 



fractional derivative is introduced on phenomenological premises |19l . |20| . In the present 
work we introduce the Levy statistics into the Langevin equation thus yielding a fractional 
FP description. This approach is similar to that of Ref. (2l| resulting in a phenomenological 



description of the non-local effects in plasma turbulence. Using fractional generalizations 



22 



m. It 



of the Liouville equation, kinetic descriptions have been developed previously 
has been shown that the chaotic dynamics can be described by using the FP equation with 
coordinate fractional derivatives as a possible tool for the description of anomalous diffusion 



24j |. Much work has been devoted on investigation of the Langevin equation with Levy white 



noise, see References 
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or related fractional FP equation [25|]. Furthermore, fractional 
derivatives have been introduced into the FP framework in a similar manner as the present 
work 
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30| but a study including drift waves is still called for. To this end we quantify 



the effects of the fractional derivative in the FP equation in terms of a modified dispersion 
relation for density gradient driven linear plasma drift waves where we have considered a case 
with constant external magnetic field and a shear-less slab geometry. In order to calculate 
an equilibrium PDF we use a model based on the motion of a charged Levy particle in a 
constant external magnetic field obeying non-Gaussian, Levy statistics. This assumption 
is the natural generalization of the classical example of the motion of a charged Brownian 
particle with the usual Gaussian statistics [31]. The fractional derivative is represented with 
the Fourier transform containing a fractional exponent. We find a relation for the deviation 
from Maxwellian distribution described by e through the quasi-neutrality condition and the 
characteristics of the plasma drift wave are fundamentally changed, i.e. the values of the 
growth-rate 7 and real frequency u are significantly altered. A deviation from the Maxwellian 
distribution function alters the dispersion relation for the density gradient drift waves such 
that the growth rates are substantially increased and thereby may cause enhanced levels of 
transport. 



The paper is organized as follows. In Sec. 2 the mathematical framework of the fractional 
FP equation (FFPE) is introduced. In Sec. 3 a dispersion relation for the density gradient 
driven drift modes using the FFPE are derived. In Sec 4 the deviations from a Maxwellian 
distribution function are investigated and the dispersion relation is solved in Sec. 5. We 
conclude the paper with a results and discussion in Sec. 6. 
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FRACTIONAL FOKKER-PLANCK EQUATION 



Following the theory of Brownian motion we write an equation of motion for a colloidal 



particle in a background medium as a Langevin equation of the following form 



3l| 



- = -.V + A(i) (1) 

Here, we assumed that the influence of the background medium can be split into a dynamical 
friction, — i/v, and a fluctuating part, A{t) which is a Gaussian white noise. The Gaussian 
white noise assumption is usually imposed in order to obtain a Maxwellian velocity distri- 
bution describing the equilibrium of the Brownian particle. This connection is due to the 
relation between the Gaussian central limit theorem and classical Boltzmann-Gibbs statis- 
tics [32]. However, the Gaussian central limit theorem is not unique and a generalization 
of the Gaussian central limit theorem to the case of summation of independent identically 
distributed random variables described by long tailed distributions is performed by Levy 



331], and Khintchine 32|. In this case the Levy distributions replace the Gaussian in a 



generalized central limit theorem. 



34| 



The simplest case of generalized Brownian motion considered by West and Seshadri 
is to assume for fluctuation part, in Equation ([T]) to be a white Levy noise. Following 



the approach used by Barkai [35|] we find the Fractional Fokker-Planck Equation (FFPE) 
with fractional velocity derivatives for shear-less slab geometry in the presence of a constant 
external force as 

ot or nis av c/v c/|v|" 
where s(= e,i) represents the particle species and < a < 2. Here, the term is 
the fractional Riesz derivative. The fractional differentiation may be represented through 
singular integrals or by its Fourier transform as we will see later in Equation (jl]). Note that 
the connection to the integral representation indicates that the model is inherently non-local 
in velocity space. The diffusion coefficient, D, is related to the damping term z/, according 



to a generalized Einstein relation 35] 



r(l + a)mr^' 

Here, Tq is a generalized temperature, and taking force F to represent the Lorentz force 
(due to a constant magnetic field and a zero-averaged electric field) acting on the particles 



of species s with mass and r(l + a) is the Euler gamma function. We find the solution 
by using the Fourier representation of equation ([2]) above as 

_£ + (_k + n.iy X 6) + z/k^)^ = -D|k^r J-., (4) 

where Vtg = CsB/nisC is the Larmor frequency of species s, 6 = B/i? is the unit vector in 
the direction of magnetic field and J^s is the characteristic function 

Fs{\^X-.t)=j y"circivexp(ik-r + ik'^- v)F,(r,v;t), (5) 

where we have denoted the wave-vector by k and the corresponding wave vector for the 
velocity as k'^. We can rewrite the kinetic equation by identification of time derivatives of 
the wave vectors as 

dt ~ dt dt (9k^ rft (9k ~ ■ ^ ' 

We use the method of characteristics on the Equation (j4]) and ([6]) whereby we find that the 



characteristics are 



^ = -D|k1"J-„ (7) 

r/k" 

— = -k + l],(k'^ X 6) + ^.k^ (8) 

0. (9) 



dt 



Following the method used in Ref. 
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30| the solution corresponding to the homogenous 
and steady state system in Fourier space is 

J•.(k^t)=e-^(l''^l'^+l''il"). (10) 



F,(r, v) = C(r) / ^ ^ ^ n "^e -l' ' r (H) 



In order to find the solution in real space we compute the inverse Fourier transform of 
Equation ffTOl) 

d\i\d\i\ 

(27r)3/2 

We define a new variable T> = ^ where coefficient D is given by the expression in Equation 
(j3]). C(r) is a normalization factor which remains to be defined. Taking the inverse Fourier 
transform of the Equation (fTTl) for a = 2 we get 

F.(r,v) = ^e-(^). (12) 



The unknown normalization factor C can be determined by comparing the integrals of the 
Maxwellian distribution and our distribution. In comparison the Maxwellian distribution is 
defined as 

F^Hv v) = e-i-l+^lWU^) (13) 



where VT,s(r) = y2Ts(r)/ms is the thermal velocity of species s. By integrating the 
Maxwellian distribution over the velocity space we find the density as 

/ rfvFf (r, v) = 2n H v^dv^ H dv^^ . ""f] e-^^+^'l^^vhir) ^ ^^^^^^ ^^^^ 

whereas performing the same integration of the expression in Equation( fT2l) we obtain 

c/vF,(r, v) = 2tt / v±dv± / dvu^-^e'^^^^ = 2n^^'^V2^C{r). (15) 

JO J-oo P 

We can now compare the two results obtained in Equations IHM - (|T5l) and we find the 
following relation 



C(r) = (16) 

The distribution function can now be determined by replacing this expression into Equation 
^ for C(r) yielding 

F,(r, v) = I^ii^) e-(^). (17) 

We can easily recover the Maxwellian distribution in Equation f lT3|) by setting a = 2 in the 
definition for D in Equation (|3]) and using that r(3) = 2. Note that for a general a, the 
equilibrium distribution is as follows 

F.(r, V) = ^^^W f ^l4^^-.(kiv,+k^v„)^-2(|.. 

27r3/2V2^i (27r)3/2 ' ^ ' 

^ = ^^77^' (19) 



where 



r(i + a 

and we have introduced a generalized thermal velocity as 



V^, = (20) 



The generalized equilibrium distribution including the effects of the fractional velocity 
derivative in Equation ( lT8l) becomes 

F.(r, v) = — ^ ^ / ^!4^^-.{k'^v,+k|fv„)^-^(|kir+|k^r)_ ^21) 



27r3/2(r(l + a))-i/2,y^y (27r)3/2 
We will now determine the dispersion relation for density gradient driven drift waves includ- 
ing the effects of the fractional velocity differential operator. 



THE DISPERSION RELATION 



In order to quantify the non-local effects on drift waves induced by the fractional dif- 
ferential operator we will determine the dispersion relation for density gradient driven drift 
modes. We start by formulating the linearized gyro-kinetic theory where the particle distri- 



bution function, averaged over gyro-phase is of the form (see Ref. [36|) 



/,(r,v) = F,(r,v) + (27r)-^x J | dk rfa;exp(ik ■ r - Jv). (22) 

We assume that the turbulence is purely electrostatic and neglect magnetic field fluctuations 
(5B = 0). For small deviations from the local equilibrium we find the linearized gyro-kinetic 
equation of the form 

(-W + k\\v\\)6f^^^{v\\,v±) + {uj- w*s)^^o(|^sr^A;±f±)i^s(a;, v)(50k,a; = 0, (23) 



where u^s = ■ is the drift wave frequency of species s, and we assumed that 

the space dependence of Fg is only in the x direction perpendicular to the magnetic field as 
well as for the density gradient. In the equation above, Jo is the Bessel function of order 
zero, f|| is the parallel velocity, v± = {v^ + VyY^'^ is the perpendicular velocity and hence we 
write the total speed as f = {v'j_ + f ^)^''^. Inserting the expression for from the Equation 
( HT^ in Equation ( 12^ and rearranging the terms we find the perturbed distribution 5/k,aj as 



Ts'kp\\-u' ' 27r3/2(r(l + a))-i/2 /2l^ 



(27r)3/2 

Here, the wave vector perpendicular to magnetic field is k± = (k"^ + kyY^"^. The gyro- 
kinetic Equation is complemented with Poisson equation for the electric potential. 
For fluctuations with wave vectors much smaller than the Debye wave vector, the Poisson 
equation becomes the quasi-neutrality condition 

Ee.5<. = 0, (25) 

s 

where the density fluctuation is related to the distribution function through 

K,. = -|^^.'^0k,. + / d^rM\n,\-'k^v^)6flJv\\,v^). (26) 

In the above equation we have separated the adiabatic response (first term on the right hand 
side) from the non-adiabatic response (second term on the right hand side). We have to keep 



in mind that the density Ug coming from the Fs{x^ v) in the adiabatic response is also given 
by Equation fl2T|) and for a general < a < 2 the adiabatic response can be different than 
that calculated by Maxwellian distribution of Equation (I13p . Using the quasi- neutrality 
condition fl2^ we find the dispersion equation which determines the eigenfrequencies as a 
function of the wave vector, uj = w(k) = ^^(k) + i7(k). In the simplest case we consider 
a plasma consisting of electrons and a single species of singly charged ions with the equal 
temperatures. For the density fluctuation therefore we have 

5nl^ = -nsir)^S<f>^AM'"''' + M^J. (27) 
^ i — I 

Therefore, the dispersion equation as in the Ref. [36( is 

^ad,e ^ ^ _M'"^'^ - M^^, (28) 

where 

J 27r3/2(r(l + a))-i/2 (27r)3/2 ^ ^ 



gives the adiabatic contribution, and 



K,a. = / ^V[ ,^ ]MbsVjVTs) X 

J rC\\V\\ CO 



1 r dkldk 



g-i(kivx+k|v||)^-P^^(ikir+|k|fr)^ (30) 



27r3/2(r(l + «))-i/2y^^7 (27r)3/2 

gives the non-adiabatic contribution. Here, hs = k±VT,s/^s- If we take a = 2 in the 
c|uation (l28!) we recover the dispersion equation for a Maxwellian distribution as in the Ref. 



36| 



Adiabatic response 

First, we may analyze the contribution from the adiabatic parts of the dispersion relation 
only by ignoring all fluctuations, yielding 

In addition, utilizing the quasi-neutrality condition while neglecting the density gradient 
in the system we have rii = n^, therefore and becomes connected through Equation 
(PTj) . This indicates that the deviation from a Maxwellian distribution described by a for 
electrons and ions becomes dependent on each other. We will get back to this relation in 
later sections. 
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DEVIATIONS FROM A MAXWELLIAN DISTRIBUTION FUNCTION 



We will now turn our attention to the problem of solving the dispersion relation described 
by E qua tion (1251) . In order to solve this dispersion equation we use the method proposed in 
361] with the difference that here we have to perform additional integrations over k''. 



Ref. 

We have 



= ^^TT^ZiQm), (32) 



where the plasma dispersion function is 

Z{Q = ^Ltm^^o r du[^p^], (33) 
with u = v\\/Vts, C,s = i^/(|^|||^s) and the function ^{v\\) is 

1 r file'" V? 

^{v\0 = ^ f ^^e-^'^^iie-TO^d'^D. (34) 



l2{r{l + a))-y^^2VS, 
The integral over vj_ can be written in a general way as 

T{bs) = 2VS, / dww^sM^v^), (35) 
where w = v^/Vts, "^s = Jo{bsV±/VTs) and, 

$(^^) = ^ / ^e-^^>-e-nfe^(l^" I"). (36) 



2(r(l + a))-y\2V^, 



(2 

The analytical solutions for integrals over with an arbitrary a in the Equations f l34|) and 
( 136|) requires rather tedious calculations. Instead we consider an infinitesimal deviation of 
the form a = 2 — e, where < e ^ 2 and expand the terms depending on a in the Equations 
( l34l) and ( l36l) around e = as follows 

e-WT^(\k I ) _ + eA(r) + 0[e% (37) 

y(r(i + a))-v2 0;^ ^Vt,s 

where 

A(r) = ^{-3 + 2-fE - 4\/t J^l' + 27E^TMfc"r 

211/4 yv^' 

+2\og[VT,s] + 2VS/og[VT,s]\k''\' + 2VSjk^\'log[\k'^\]}. (38) 



Here, we have used the Euler-Mascheroni constant 7^; = 0.57721. The first term in Equation 
f l37|) will produce 



*W = :fT372' and $(«;)= (39) 

which give the Maxwellian adiabatic response 

W"^'' = 1. (40) 

By using the expansion defined by the expression fISTI) in Equations fl2^ and fl2^ . the 
adiabatic and non-adiabatic part of the dispersion relation M°''^'^ and ^ are as follows 



/OO /"OO 
-00 JO 



2y27r3/2 7 (27r)3/2 



^ II e-^^*^" ^-+'^i^")A(fcl)A(A:|))6 + = 1 + eW^"^' . 



(41) 



and 



/OO i-OO /J /j! 

-00 Jo K\\V\\ — UJ 



2V27r3/2 7 (27r)3/2 



(v^^T,.(r))3 

/OO (-00 /J — 
dvii / cit;xt^±[^ ^]^s(&st^±/^Ts) X 
-00 Jo k\\v\\ — UJ 

^^e-^(^^^-+'^^ll)A(A:l)A(fc|))e + 0[ef = N^^ + eW^^ 



(42) 



Inserting these relations we may rewrite the dispersion relation f l28|) in the form 

(1 + NIJ + eiW^"^^ + WIJ = -(1 + NIJ - e{W^'''^ + WIJ. (43) 
The first terms on the right and left hand sides generate the usual contributions to the dis- 



36| and the terms proportional to e generate the non-Maxwellian 



persion equation as in Ref. 
contributions where we have 

= ^^TlV^ZiQm), (44) 
\I^\\\Vt,s 

with the usual plasma dispersion function Z{^s) written as 

Zi^s) = ^Ltm^^o / due-'' [ -], (45) 

VTT J-00 M — 4s — "^CT 

10 



and 

T{hs) = 2 / dwiue-"" (46) 
Jo 

The effects of the fractional velocity derivative can be boiled down to a non-Maxwellian 
contribution of the form 

Wi^ = ^lfr^Z,{QTM, (47) 
where the non-Maxwellian plasma dispersion function is given by 

Z.i^s) = ^Lim^^o / du[ ^L], (48) 

a/tt j-oo u — t,s — icr 

with the function ^{v\\) being 

1 r dk" 

= ^.-M-^^v,OK{^). (49) 

It is important to note that the deviation from Maxwellian is different for the different 
species (electrons and ions). In the rest of Sec. 4, we will quantify the deviations. The 
non-Maxwellian contribution to Equation (135|1 is 

/■oo 

T,{hs) = 2VS, / dww^HsM^v^), (50) 

where 

<^>K) = ^ / ^exp(-zA:>^)A(A:l). (51) 

To extimate the non-Maxwellian contribution we need to determine the inverse Fourier 
transforms of the Equations (l49l) and (|5T1) resulting in 

1 



|-4(-2 + ^e)z' + (-7 + 4 7s) + 2 log[VT,.] + 2e^\Fi[-, -z"]^ (52) 

with z = {u^w} and iFi[a;6;z] denoting Kummer's confluent hypergeometric function. 
Therefore we can write 

y^^ad,s ^ / / (53) 

TT J-oo JO 



By inserting typical values for the plasma parameters from Ref. jol we find the velocities as 
Vr.e = 5.93 X 10^[cm/s] and Vt^i = 1.38 x l{f[cm/ s] and we obtain 

^ad,e ^ 33 724 ^ w'"^'' = 23.6591. (54) 
11 



Following the adiabatic condition in Equation ( 13T]) and the expanded dispersion relation in 
Equation fH3|) we obtain the following ratio between the non-Maxwellian contributions 

- = jjj-r = 1.42541. (55) 

This relation means that if there is a deviation of the distribution function from the 
Maxwellian for plasma electrons, the deviation from the Maxwellian for ions will be ~ 1.4 
larger. 



SOLUTIONS OF THE DISPERSION RELATION 

We will solve the dispersion relation in terms of expansions of the plasma dispersion 
function by noting that the drift waves are defined in the frequency range |A;|||VTj ^ w ^ 
|fc|||VTe in evaluating Equations ( 115|) and ( 118|) . We define the expansion parameter for 
electrons in powers of = I^Te) ^ 1 and for ions we expand it in powers of ^"^ = 

(I |VTi)/w ^ 1, respectively. The Maxwellian dispersion function Zi^^^ has the same 
definition as in Ref. 



Z{Q = -^Lzm^^o / due-"^ [ -] = -2^^ + ^+ ^V^{1 - C') + 

a/TT J-oo u — t^e — la 6 

whereas the non-Maxwellian plasma dispersion function Z^{C,e) becomes 

a/vt j-oo m — 4e — "^cr y/n 



du<!>{u)[^— + ^\ + , ^3 + . + 0[C]]- (57) 

3 u — la [u — lay [u — lay [u — lay 

For ions, using the expansion in powers of we can rewrite the above integrals as a 
function of the expansion parameter as 

Z(6) = ^Lim^^o / due-- [ -] = -^^ - + ^[C'], (58) 

a/vt j-oo u — t,i — icr 2 



and the non-Maxwellian Z^{C,i) becomes 



Ze(^i) = ^Lim^_>o / du[ — ^^^^ . ] = -^Lim, 



/tt 7-00 u — — ia ^vr 

1 u u"^ 



o— i>0 
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We can now evaluate he Maxwellian integrals of the forms T{be) and T{bi) assuming \&e = 1, 
^» = Jo{biV±/VTi) we get 



T{be) = 2 / dwwe-"" = 1, (60) 
Jo 

and 

r(6,) = 2 / dwwe-"" <ile{kw) = e-'^/^Ioik), (61) 
Jo 

where Xq denotes modified Bessel function of the zeroth order. The final result will be found 
after evaluating the non-Maxwellian T^{bf.) and T^{bi) are given as 

POO 

r,(6e) = 2VSe / dww<!>{w) = 4.8 x 10^ (62) 

and 

roo 

r,{bi) = 2VSi / dww'^{biw)<^{w) = 6.1 x 10^ (63) 
' Jo 

where we have used Vr,e = 5.93 10^[cm/s], Vr.i = 1-38 10^[cm/s] and bi = 0.1. Finally we 
can summarize different terms in the dispersion relation (l43l) as 



Wl^ = i^,-u.,)Z,i^,)TM, (64) 

where a)*,.) = a;*_s/|^|| | Vt,s- Note that the non-Maxwellian contributions in Equations (1571) . 
(!59|) . fl62l) and (!63l) have been calculated numerically. By utilizing the found values of the 
integrals above we rewrite the dispersion relation P3l) as follows 



(1 + e^W'^n + (Ce - C^*,e){-2^e + ^ + ^V^(l " C') + ^.Z ,Ue)T ,ib,)} = 

_(1 + e.iy'^^'O - (e. - - ^C')e-'''/'Xo(60 + e.Z,(e0r.(60} (65) 

where W"''^'^ are given in Equation (15^ and we will use the ratio between eg and from 
Equation (j55|) . 
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RESULTS AND DISCUSSION 



We have derived a dispersion relation for drift waves driven by a density gradient in 
a shear-less slab geometry with constant magnetic field where the small deviation from a 
Maxwellian distribution is described by e. Here we will determine the quantitative effects 
on the real frequency and growth rate as a function of this deviation. We start by assuming 
that we have adiabatic electrons for which the dispersion Equation fl65l) is, 

2 + 6,(2.35 w"^''' + W"^^') = 
-(6 - io^M-^i' - \^-')e~'^^'Ubi) + e,Z,{^,)TM}. (66) 

After rearranging the terms in the above equation we finally get the following relation for 

^ -2ef + [ef + 0-5^i - - 0.5cI;,,]e-^'/%(&0 

where W'^'^'^"* = 2.35 W^^'^ + W^^''. This relation gives the possible deviation of the equi- 
librium PDF from the Maxwellian PDF for a given plasma turbulence, i.e ^j. One has to 
remember that only positive values of Re[e] are physically meaningful. 

Using the same plasma parameters as was used in Equations ( l5l|) and ( |6T][63|) we compute 
the term Z^{C,i), and from Equation f l59|) we get 

Vt ■ 1 1 

Z,(e^) = -^Ltm^^oi— — / du<^{u) + — — / u^du<!>{u)} 

Vvr (— 4i — ta) J-oo (— 4j — ^c^) J-oo 

-6.5 X 10-" - 3.8 X 10-"^ 
= g ■ (68) 

Here, those integrations omitted resulted in zero contributions and rewriting Equation ( 166|) 
by using these explicit values results in the expression for the deviation in Equation fl67|) we 
obtain 

^ -2C' + i^f + 0.5^. - - 0.5a;,,je-^'/%(60 

66.3^3 - 39.2^^ + 39.2a;,,i + 23.0(x;,42 ) 

Figure [T] shows from Equation ( 1^ where = u + i'j. Here, the values of co,'y are 
normalized to |fc|||VT,i. We have assumed parameter values 6j = 0.1, fcy = 10"'^ and 0)*^^ = 
—7.1 X 10^ with d In n/dx = 1. It is found that there is a threshold in the growth rate 7 
close to 7 = 0.7 and that increasing to 1.0 only increases the deviation from a Maxwellian 
from to 0.03. It should be noted that e increases the excess kurtosis of the distribution 



14 



1.0 



FIG. 1. e as a function of uj and 7 where = cj + i^. We have assumed hi = 0.1, /cy = 10 ^ and 
Lj^^j = —7.1 X 10^ with d In n/dx = 1. 

function by a similar amount thus a quite small deviation from a Maxwellian can have a 
rather significant impact. 

In figure Ej the mode growth rate as a function of is shown. Note that in this figure the 
values of growth rate are the solutions of the Equation ( l69l) for a given while in the figure 
[T]we solve Equation fl69|) for e at a given ^j. As the dispersion equation is of 3rd order in cu 
three possible solutions exist, however we are only interested in the solutions with non-zero 
imaginary value, 7 > corresponding to unstable situations. It is shown in figure [2] that a 
deviation of = 0.01 yield an increase of about 20% in the growth rate. Furthermore, the 
growth rate increases almost linearly with increasing and such an increase in the growth 
rate will lead to a significant increase in the level of anomalous fiux. 

In summary, we have derived a dispersion relation for density gradient driven linear drift 
waves including the effects coming from the inclusion of a fractional velocity derivative in the 
Fokker-Planck equation in the case of constant magnetic field and a shear-less slab geometry. 
The solutions of this Fokker-Planck equation are the alpha-stable distributions. It has not 

alt been shown that in a direct way one can derive the alpha-stable distribution function 
, l37| from the classical form of collision operator 38|. One way may be to construct a 



new type of coUisional operator by considering a fractal phase space and reformulate the 

15 
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FIG. 2. 7 as a function of e. The same plasma parameters as in the figured] are used. 

collision operator on this new space. However, such a discussion is outside the scope of 
the present paper. Interestingly enough, we note that non-local effects are observed in 
non-linear coUisionless fluid simulations of plasma turbulence where the non-local transport 
showing Levy features are induced by the interaction of the non-linear terms in the dynamical 



equations 
analytica 
fluxes jl6 |. 



12| . The non-local features of non-linear fluid models are indicated by recent 



theories using path-integral methods to derive probability density fucntions of 



The fractional derivative is represented with the Fourier transform containing a fractional 
exponent that we are able to connect to the deviation from a Maxwellian distribution de- 
scribed by e. The characteristics of the plasma drift wave are fundamentally changed, i.e. 
the values of the growth-rate 7 and real frequency u are significantly altered. A deviation 
from the Maxwellian distribution function alters the dispersion relation for the density gra- 
dient drift waves such that the growth rates are substantially increased and thereby may 
cause enhanced levels of transport. 
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